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ABSTRACT 


In Part I a method for calculating wave packet 
trajectories and wave heights is based on the assumption 
that the water depth contours are locally parallel in the 
vicinity of each ray point. This method is extended in 
order to predict the modification to surface gravity water 
waves in shoaling water when the water depth contours are 
not parallel. The calculations are greatly simplified by 
choosing a coordinate system at each ray point in which one 
axis is aligned parallel with the direction of the gradient 
of the water depth. Example printouts and plots are 
presented to illustrate the wave prediction method. It is 
discovered that when waves initially approach sinuous water 
depth contours symmetrically with respect to the beach 
there can be more energy in the bays than at the headlands. 
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CHAPTER I INTRODUCTION 


This is a supplement to the report by Breeding, 
J. Ernest, Jr., K.C. Matson, and Nourollah Riahi, "A Method 
for Calculating Wave Packet Trajectories and Wave Heights," 
Department of Oceanography, Florida State University, 
Tallahassee, (Report No. JEB-1 (1978) hereafter referred to 
as_Pant—l,, e program for predicting the modification to 
“waves in shoaling water described in Part I is developed on 
the assumption that the water depth contours are locally 
parallel in the vicinity of each ray point. This method 
results in accurate wave predictions for many examples of 
bottom topography. However, the wave height predictions are 
not accurate when the water depth contours are highly non- 
parallel in the region of a ray point. 


The objective in this report is to extend Part I in 
order to make accurate wave predictions when the water depth 
contours are not parallel. The theory is developed in 
Chapter II. The calculations are greatly simplified by 
making them in a x'y'-coordinate system which is chosen so 
that at each ray point the positive x'-axis is in the direc- 
tion of the gradient of the water depth. As a result, the 
first partial derivatives of the water depth and wave 
velocities with respect to y' vanish. Further, there is a 
simplification in the second partial derivatives involving 
y'. In Chapter III modifications to the computer program 
are described. To illustrate the wave prediction program 
two examples of nonparallel water depth contours are pre- 
sented in Chapter IV. An interesting result is that when 
waves initially approach sinuous water depth contours 
symmetrically with respect to the beach (headon) it is 
possible for there to be more energy in the bays than at 
the headlands. This is in contrast to the result for mono- 
chromatic trajectories where under similar conditions there 
is always more energy predicted at the headlands than in 
the bays. 


CHAPTER II ‘THEORY FOR NONPARALLEL WATER DEPTH CONTOURS 


2.1 Spatial Derivatives of the Water Depth. At each 
point of a wave packet trajectory the calculations are made 
in a x'y'-coordinate system where the x'-axis is taken in 
the direction of the gradient of the water depth. The par- 
ticulars of a wave packet trajectory are tabulated ina 
xy-coordinate system which retains a fixed orientation with 
respect to the water depth grid. The relationships between 
these coordinate systems and a specific ray point for a set 
of nonparallel water depth contours are shown in Figure (2-1). 
Equations relating these coordinate systems are given by 
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where a is the angle by which the x'-axis is rotated with 
respect to the x-axis and h is the water depth. 


The partial derivatives of h in the x'y'-coordinate 
system with respect to the partial derivatives of h in the 
xy-coordinate system are given by 
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Figure (2-1). Relationships between the coordinate 
systems and the water depth contours. 
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In Equation (2-5), dh/day' = 0 as a result of choosing the 
positive x'-axis in the direction of the gradient of the 
water depth. If the water depth contours are locally paral- 
lel, which was the case dealt with in Part I, then 
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2.2 Spatial Derivatives of the Wave Speeds. In this 
section relations are presented for connecting the geometric 
group speed G, the collinear group speed U, the phase speed 
v, and h. The water depth contours are assumed to be non- 
parallel. 


a. Derivatives of v 


The phase speed of a surface gravity water wave is 
defined (Part I) 
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where T is the wave period and g is the acceleration due to 
gravity. 


The first partial derivatives of v in the x'y'-coordinate 
system are given by 
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The second partial derivatives of v are defined by 
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b. Derivatives of U 


The collinear group speed of a surface gravity water 
wave is defined (Part I) 
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The first partial derivatives of U in the x'y'-coordinate 
system are given by 
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The second partial derivatives of U are given by 
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ec. Derivatives of G 
The geometric group speed is defined by 
= U wm (2-33) 


where 


(2-34) 


in which 6' is the direction of the wave packet, y' is the 
direction of the wavelets, and both directions are measured 
with respect to the positive x'-axis. The first spatial 
derivatives of G are 
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The ray curvature for a wave packet is given by (Part I) 
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where ds, is an element of arc length along the packet tra- 
jectory. Since the first derivatives in y' vanish it is found 
that 
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In like manner, when considering the ray curvature of the 
wavelets it is found that 
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When Equation (2-42): is substituted into (2-35) and the 
result is pag it is found that 
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The second partial derivative of G with respect to x' is 
given by 
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The remaining second order spatial derivatives of G are 
obtained by differentiating Equation (2-36). Since the first 
derivatives in y' vanish it is found that 


6 _ FU oor 2° 
@y) = Dp Od U dm oy (2-50) 


< 8e Cord - U Aun > 


(2-51) 
ao ay! ax'ay! 


The second partial derivative of » with respect to y' is 
found by differentiating Equation (2-43). The result is 


(2-52) 


When Equation (2-52) is substituted into (2-50) and the terms 
rearranged it is found that 
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After differentiating Equation (2-43) with respect to x' it 
is found that 
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When Equation (2-56) is substituted into (2-51) the simpli- 
fied result becomes 


(2-57) 


2.3 Ray Curvature for Nonparallel Water Depth Contours. 
The ray curvature kK, depends only on derivatives of the first 
order. By making the calculations in the x'y'-coordinate 
system the first derivatives in y' vanish. This results in 
a simplified ray curvature expression which is formally the 
same as presented in Part I for locally parallel water depth 
contours. 
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2.4 Ray Separation Equation for Nonparallel Water Depth 
Contours. The ray separation equation is defined (Part I) 


44 + et + 46 =0 (2-59) 


where 8 is the ray separation factor, t is time, and 
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In the x'y'~coordinate system p is simplified as is Kq In 
q there is also a reduction in the cre for the 
second order derivatives involving y' 


CHAPTER III SUBROUTINE SURFCE 


3.1 Modifications to Subroutine SURFCE. Modifications 
have been made to the program subroutine SURFCE in Part I in 
order to remove the restriction that the water depth contours 
be locally parallel about each ray point and to simplify some 
of the computations. When the program is run with these 
modifications minor changes will occur in the printed output 
for the sample input data presented in Part I. 


Referring to the program listing for SURFCE on pages 
52-54 of PartI the following modifications have been made. 
In place of SURFACE 58 an alternative expression is used to 
compute dh/dx' (DHDX). Between SURFACE 87 and SURFACE 88 
dI/dx' (DIDX) is computed as are the hyperbolic sine, cosine, 
and tangent of I. Statements SURFACE 91 through SURFACE 93 
for computing 3U/ax' (DUDX) have been replaced by a simpli- 
fied expression. Between SURFACE 107 and SURFACE 108 
a7h/(dy')? (DHDYY) and 3*h/ax'dy' (DHDXY) are computed. 
SURFACE 111 has been replaced by expressions to compute 
d*v/(dy')* (DVDYY) and 3*v/dx'dy' (DVDXY). Statements 
SURFACE 114 through SURFACE 120 have been replaced by expressions 
for calculating 07I/(dy')? (DIDYY) and 37I/dx'dy' (DIDXY). A 
simplified expression has been substituted for 37U/(dx')? (DUDXX). 
Also, the values of 97U/(ay')* (DUDYY) and 97U/dx'dy' (DUDXY) 
are computed. Between SURFACE 122 and SURFACE 123 ¢ (ZETA), 
& (XI), 37G/(day')? (DGDYY), and 37G/dx'dy' (DGDXY) are 
determined. Simplified expressions are used to compute 
a?G/(day')? and 9?7G/ax'dy' when!|tan 6'| < tan 5° and |tan y'|< tan 5°. ' 
Statement SURFACE 124 has been replaced by the complete 
expression for computing q (QOT). 


3.2 Description of Subroutine SURFCE. SURFCE is called 
by RAYN and MOVE to calculate h, a, y, G, p, q, kg, and other 
ray particulars. At the first ray point twelve values of h 
from CMAT are selected about the point as shown in Figure (2-2) 
of Part I. A quadratic surface is fit to the set of water 
depths. At successive ray points the quadratic surface is 
determined only if there is a change in the set of twelve 
water depths. The water depth and its partial derivatives in 
the fixed xy-system, 3h/3x, dh/dy, 37h/ax*, 37h/dy*, and 
d*h/dxdy, are determined at the ray point by interpolating on 
the quadratic surface. 


If h < 0, NDP = 2 and there is a RETURN. If h > 0O the 
ratio of the water depth to the deep water wavelength is 
computed. If h/d\g > 0.64, which defines deep water, NFK = l. 
If h/Aq s 0.64, NFK = 2, VELCTY is called, and after the 
return if NFK = 1, W= 0. If NFK = 2, CONDER is called to 
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13 ? 


compute W. The values of dv/dx and dv/dy are calculated 
using W. 


At each ray point calculations are made ina x'y'- 
coordinate system which is chosen with the positive x'-axis 
in the direction of the gradient of the water depth. The 
value of dh/dx' is computed, and if it exceeds 0.00001 the 
angle a by which the x'-axis is rotated with respect to the 
x-axis is computed. If | dah/ax'| < 0.00001 the water depth is 
assumed to be constant and a remains constant. 


If FLAG1 = 0, y' is computed, and if necessary it is 
placed within the range |y'| < 360°. A check is made to 
determine if there is total reflection. If there is, 

FLAG2 = 1 and there is a RETURN. Otherwise, FLAG2 = 0 and the 
new y' is computed using Snell's law with phase velocity 
following a set of rules. Using the values of y', y is 
computed. When FLAG 1 # 0 these steps for computing the new 
wavelet direction and the test for total reflection are 
omitted. 


The values of »¢ , G, and dv/dx' are calculated. If 
NFK = 2, dU/dx' is determined using its unsimplified expres- 
sion. If NFK # 2, the deep water formula is used to calculate 
dU/ax'. The value of 3U/3x' is used in computing 3G/ax'. 


If NFK # 2, the coefficients of the ray separation 
equation and the ray curvature are set equal to zero. Then 
there is a RETURN. If NFK = 2, p, d7*h/(dx')?, a*h/(ay')®, 
d*h/dx'dy', d7v/(dx')?, a2v/(dy' o*v/ax'dy', 37U/(dx')*, 


d7U/(dy')*, d7U/dx'dy', 37G/(ax')?, 37G/(day')?, a2G/(axtay'), 
q, and K, are computed. This is followed by a return. 
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CHAPTER IV WAVE PREDICTION EXAMPLES FOR 


NONPARALLEL WATER DEPTH CONTOURS 


Two examples of nonparallel water depth contours are 
used to illustrate the wave prediction program. 


4.1 Sinusoidal to Parallel Water Depth Contours. 
Figure (4-1) shows a set of wave packet trajectories for a 
wave period of 7.0 seconds. The rays begin in deep water 
(initial water depth equal to or greater than one half the 
wavelength). The water depth contours are sinusoidal at the 
shoreline and gradually become parallel at a water depth of 
38.6 meters. The amplitude of the water depth contour at the 
shoreline is 2.5 kilometers, and the contour wavelength is 
10 kilometers. GRID, the distance between grid points, is 
156.25 meters. 


From the figure it is seen that the energy is fairly 
evenly distributed along the coastline with slightly more 
energy at the headland than in the bay. Figures (4-2) through 
(4-4) show the printed output for rays number 4, 8, and 16, 
respectively, of Figure (4-1). The computed refraction 
coefficients are in good agreement with values estimated from 
the plot. 


4.2 Sinusoidal Water Depth Contours. In Figure (4-5) 
the period of the wave packet trajectories is 14.0 seconds 
and they begin in deep water. The water depth contours 
are sinusoidal with an amplitude of 5 kilometers and a 
wavelength of 20 kilometers. GRID has a value of 312.5 
meters. 


This example is quite interesting since there is decidedly 
more energy in the bay than at the headland. The opposite 
result would be expected for monochromatic trajectories. The 
refraction of wave packets could explain why there is more 
erosion in bays than at headlands for some coastlines. 

Figures (4-6) through (4-8) show the printed output for rays 
number 4, 8, and 16, respectively, of Figure (4-5). The 
computed refraction coefficients agree favorably with values 
estimated from the plot. 
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Figure (4-1). Plot for sinusoidal to parallel water 


depth contours. 
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Figure (4-5). Plot for sinusoidal water depth contours. 
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Figure (4-5). 
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